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Abstract 

The micromaser possesses a variety of dynamical phase transitions parametrized by the 
flux of atoms and the time-of-flight of the atom within the cavity. We discuss how these 
phases may be revealed to an observer outside the cavity using the long-time correlation 
length in the atomic beam. Some of the phase transitions are not reflected in the average 
excitation level of the outgoing atom, which is the commonly used observable. The corre- 
lation length is directly related to the leading eigenvalue of the time evolution operator, 
which we study in order to elucidate the phase structure. We find that as a function 
of the time-of-flight the transition from the thermal to the maser phase is characterized 
by a sharp peak in the correlation length. For longer times-of-flight there is a transition 
to a phase where the correlation length grows exponentially with the flux. We present a 
detailed numerical and analytical treatment of the different phases and discuss the physics 
behind them. 
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1 Introduction 



The highly ideahzed physical system of a single two-level atom in a superconduct- 
ing cavity, interacting with a quantized single-mode electromagnetic field, has been 
experimentally realized in the micromaser and microlaser systems 0. De- 

tails and references to the literature can be found in e.g. the reviews 0-||l3l. In 
the absence of dissipation (and in the rotating wave approximation) the two-level 
atom and its interaction with the radiation field is well described by the Jaynes- 
Cummings (JC) Hamiltonian [14]. Since this model is exactly solvable it has played 
an important role in the development of modern quantum optics (for a recent ac- 



count see e.g. Refs. []T2|, [T^). The JC model predicts non-classical phenomena, such 
as revivals of the initial excited state of the atom |jl5|-||2^, experimental signs of 
which have been reported . 



Correlation phenomena are important ingredients in the experimental and theo- 
retical investigation of physical systems. Intensity correlations of light was e.g. used 
by Hanbury-Brown and Twiss p2[ as a tool to determine the angular diameter of 
distant stars. The quantum theory of intensity correlations of light was later devel- 



oped by Glauber These methods have a wide range of physical applications 

including investigation of the space-time evolution of high-energy particle and nuclei 
interactions P3, ES. In the case of the micromaser we have recently suggested 



that correlation measurements on atoms leaving the micromaser system can be used 
to infer properties of the quantum state of the radiation field in the cavity. 

In this paper we present a detailed account of the role of long-time correlations in 
the outgoing atomic beam and their relation to the various phases of the micromaser 
system. Fluctuations in the number of atoms in the lower maser level for a fixed 
transit time r is known to be related to the photon-number statistics |p7|-||30 



The experimental results of |^ are clearly consistent with the appearance of non- 
classical, sub-Poissonian statistics of the radiation field, and exhibit the intricate 
correlation between the atomic beam and the quantum state of the cavity. Related 
work on characteristic statistical properties of the beam of atoms emerging from the 
micromaser cavity may be found in Ref. 

The paper is organized as follows. In Section || we discuss the standard theo- 
retical framework for the micromaser and introduce some new notation. A general 
discussion of long-time correlations is given in Section ^ where we also determine 
the correlation length numerically. Before entering the analytic investigation of the 
phase structure we introduce some useful concepts in Section ^ and discuss the eigen- 
value problem for the correlation length. The heart of the paper lies in Section |^, 
where details of the different phases are analysed. In Section ^ we study effects re- 
lated to the finite spread in atomic velocities. The phase boundaries are defined in 
the limit of an infinite fiux of atoms, but there are several interesting effects related 
to finite fiuxes as well. We discuss these issues in Section 0. Finally we summarize 
our results in Section R|. 
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2 Basic micromaser theory 



In the micromaser a beam of excited atoms is sent through a cavity and each atom 
interacts with the cavity during a well-defined transit time r. The theory of the 



micromaser has been developed in |27, 28 1, and in this section we briefiy review the 



standard theory, generally following the notation of that paper. We assume that 
excited atoms are injected into the cavity at an average rate R and that the typical 
decay rate for photons in the cavity is 7. The number of atoms passing the cavity in 
a single decay time N = R/'j is an important dimensionless parameter, effectively 
controlling the average number of photons stored in a high-quality cavity. We shall 
assume that the time r during which the atom interacts with the cavity is so small 
that effectively only one atom is found in the cavity at any time, i.e. Rt <^ 1. A 
further simplification is introduced by assuming that the cavity decay time I/7 is 
much longer than the interaction time, i.e. 7r -C 1, so that damping effects may be 
ignored while the atom passes through the cavity. This point is further elucidated 
in Appendix ^ In the typical experiment of Ref. these quantities are given the 
values = 10, Rr = 0.0025 and 7r = 0.00025. 



2.1 The Jaynes— Cummings Model 

The electromagnetic interaction between a two-level atom with level separation ujq 
and a single mode with frequency u) of the radiation field in a cavity is described, in 



the rotating wave approximation, by the Jaynes-Cummings (JC) Hamiltonian |T4| 



H = uja*a + ^ijJqOz + g{aa+ + a*cr_) , (2.1) 



where the coupling constant g is proportional to the dipole matrix element of the 
atomic transition^. We use the Pauli matrices to describe the two-level atom and 
the notation a± = {(Tx ± iay)/2. For g = the atom-plus-field states \n,s) are 
characterized by the quantum number n = 0, 1, ... of the oscillator and s = ± for 
the atomic levels (with — denoting the ground state). At resonance u = ujq the 
levels \n — 1, +) and \n, — ) are degenerate for n > 1 (excepting the ground state 
n = 0), but this degeneracy is lifted by the interaction. For arbitrary coupling g 
and detuning parameter Aa; = uq — u the system reduces to a 2 x 2 eigenvalue 
problem, which may be trivially solved. The result is that two new levels are formed 
as superpositions of the previously degenerate ones with a separation in energy 
~ -E'n,- = y/'^uj'^ + 4:g'^n. The system performs Rabi oscillations with this 

-•^This coupling constant turns out to be identical to the single photon Rabi frequency for the 
case of vanishing detuning, i.e. g ~ fl. There is actually some confusion in the literature about 
what is called the Rabi frequency |Q. With our definition, the energy separation between the 
shifted states at resonance is 2Q. 
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frequency between the original, unperturbed states with transition probabihties 



(n-l,+|e-^^^|n,-)p = g„(r) , 
(n,+|e-^^^|n,+)|2 = 1 - g„+i(r) 
(n + l,-|e-^^^|n,+)|2 = g„+i(r) . 



(2.2) 



These are all expressed in terms of 



Notice that for Acj = we have = sin^((7rA/n). Most of the following discussion 
will be limited to this case. 

Denoting the probability of finding n photons in the cavity by pn we find the 
conditional probability that an excited atom decays to the ground state in the cavity 
to be 

oo 

^(-) = (Wl) = E WlPn . (2.4) 
n=0 

It is this sum over the incommensurable frequencies g^/n that is the cause of some 
of the most important properties of the micromaser, such as quantum collapse and 
revivals (see e.g. Refs. 01-|3^). These effects are most easily displayed in the case 



that the cavity field is coherent with Poisson distribution 

Pn = ^e-<") . (2.5) 

In the more realistic case, where the changes of the cavity field due to the passing 



atoms is taken into account, a complicated statistical state of the cavity arises |27 
3^-[B3] It is the details of this state that are investigated in this paper. 



2.2 Mixed states 

The above formalism is directly applicable when the atom and the radiation field 
are both in pure states initially. In general the statistical state of the system is 
described by an initial density matrix p, which evolves according to the usual rule 
p — i> pit) = exp{—iHt)pexp{iHt). If we disregard, for the moment, the decay of 
the cavity field due to interactions with the environment, the evolution is governed 
by the JC Hamiltonian in Eq. (171). It is natural to assume that the atom and the 
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radiation field of tlie cavity initially are completely uncorrelated so that the initial 
density matrix factorizes in a cavity part and a product of k atoms as 



P = PC ® PAi ® PA2 ® ■ ■ • ® PAfc ■ (2.6) 

When the first atom Ai has passed through the cavity, part of this factorizability is 
destroyed by the interaction and the state has become 

p{r) = pcAi (^) ® PA2 ® ■ ■ ■ ® PAfe • (2.7) 

The explicit form of the cavity-plus-atom entangled state Pc,Ai(t) is analysed in 
Appendix |^. After the interaction, the cavity decays, more atoms pass through and 
the state becomes more and more entangled. If we decide never to measure the state 
of atoms Ai . . . Ai with i < k, we should calculate the trace over the corresponding 
states and only the po-component remains. Since the time evolution is linear, each 
of the components in Eq. ( p. 71) evolves independently, and it does not matter when 
we calculate the trace. We can do it after each atom has passed the cavity, or 
at the end of the experiment. For this we do not even have to assume that the 
atoms are non-interacting after they leave the cavity, even though this simplifies the 
time evolution. If we do perform a measurement of the state of an intermediate 
atom Ai, a correlation can be observed between that result and a measurement of 
atom Ak, but the statistics of the unconditional measurement of A^ is not affected 
by a measurement of Ai. In a real experiment also the efficiency of the measuring 
apparatus should be taken into account when using the measured results from atoms 
Ai, . . . ,Ai to predict the probability of the outcome of a measurement of Ak (see 



Ref. ||3^ for a detailed investigation of this case). 

As a generic case let us assume that the initial state of the atom is a diagonal 
mixture of excited and unexcited states 




where, of course, a,b > and a+b = 1. Using that both preparation and observation 
are diagonal in the atomic states, it may now be seen from the transition elements in 
Eq. (|2.2|) that the time evolution of the cavity density matrix does not mix different 
diagonals of this matrix. Each diagonal so to speak "lives its own life" with respect 
to dynamics. This implies that if the initial cavity density matrix is diagonal, i.e. of 
the form 

oo 

PC = J2P'^\^)i^\ ' (2-9) 

n=0 
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with Pn > and J2'^=oPn = 1; then it stays diagonal during the interaction between 
atom and cavity and may always be described by a probability distribution Pn{t)- 
In fact, we easily find that after the interaction we have 



p„(r) = ag„(r)p„_i + 6g„4.i(r)p„+i + (1 - ag„+i(r) - 6g„(r))p„ , (2.10) 



where the first term is the probability of decay for the excited atomic state, the 
second the probability of excitation for the atomic ground state, and the third is 
the probability that the atom is left unchanged by the interaction. It is convenient 



to write this in matrix form 55 



p(r) = M(r)p , (2.11) 

with a transition matrix M = M(+) + M(— ) composed of two parts, representing 
that the outgoing atom is either in the excited state (+) or in the ground state (— ). 
Explicitly we have 

M{ + )n.m = bqn+lSn+l,m + «(! " qn+l)Sn,ni , 

(2.12) 

M{-)nm = aqnSn,m+l + " gn)^n,m • 

Notice that these formulas are completely classical and may be simulated with a 
standard Markov process. The statistical properties are not quantum mechanical 
as long as the incoming atoms have a diagonal density matrix and we only measure 
elements in the diagonal. The only quantum mechanical feature at this stage is 
the discreteness of the photon states, which has important consequences for the 
correlation length (see Section p.3|) . If the atomic density matrix has off-diagonal 
elements, the above formalism breaks down. The reduced cavity density matrix will 
then also develop off-diagonal elements, even if initially it is diagonal. We shall not 
go further into this question here (see for example Refs. ||40|| - 



2.3 The lossless cavity 

The above discrete master equation ( |2.10| ) describes the pumping of a lossless cavity 
with a beam of atoms. After k atoms have passed through the cavity, its state has 
become M'^p. In order to see whether this process may reach statistical equilibrium 
for oo we write Eq. ( |2.10| ) in the form 



Pn{r) =Pn + Jn+1 - Jn , (2.13) 
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where J„ = —aqnPn-i + bqnPn- In statistical equilibrium we must have J„+i = J^, 
and the common value J = Jn for all n can only be zero since p„, and therefore J, has 
to vanish for n — > oo. It follows that this can only be the case for a < b i.e. a < 0.5. 
There must thus be fewer than 50% excited atoms in the beam, otherwise the lossless 
cavity blows up. If a < 0.5, the cavity will reach an equilibrium distribution of the 
form of a thermal distribution for an oscillator pn = {1 — a/b){a/b)^. The statistical 
equilibrium may be shown to be stable, i.e. that all non-trivial eigenvalues of the 
matrix M are real and smaller than 1. 



2.4 The dissipative cavity 

A single oscillator interacting with an environment having a huge number of degrees 
of freedom, for example a heat bath, dissipates energy according to the well-known 
damping formula (see for example [0, 

i[pc,uja*a\ 

-|7(^6 + l){a*apc + pca*a - 2apca*) (2-14) 
-^'ynb{aa*pc + pcaa* - 2a*pca) , 

where rib is the average environment occupation number at the oscillator frequency 
and 7 is the decay constant. This evolution also conserves diagonality, so we have 
for any diagonal cavity state: 




= -inb + l)inpn- in + l)pn+i) -nh{{n + l)pn-npn-^i) , (2.15) 

7 dt 

which of course conserves probability. The right-hand side may as for Eq. (|2.13|) be 
written as Jn+i — Jn with J„ = {rib + l)^Pn — nbupn-i and the same arguments as 
above lead to a thermal equilibrium distribution with 

p„ = ^f^V . (2.16) 



1 + nfc Kl + rib 



2.5 The discrete master equation 



We now take into account both pumping and damping. Let the next atom arrive in 
the cavity after a time T ^ r. During this interval the cavity damping is described 



by Eq. (|2.15|) , which we shall write in the form 



dp 

- = -^Lcp 



(2.17) 
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where Lc is the cavity decay matrix from above 



{Lc)nm = {nb+l){nSn,rn- {n+l)Sn+l,m) +nb{in+l)Sn,m-nSn-l,m) ■ (2.18) 

The statistical state of the cavity when the next atom arrives is thus given by 

p{T) = e-^^''^M{T)p . (2.19) 

In using the full interval T and not T — r we allow for the decay of the cavity in 
the interaction time, although this decay is not properly included with the atomic 
interaction (for a correct treatment see Appendix 

This would be the master equation describing the evolution of the cavity if the 
atoms in the beam arrived with definite and known intervals. More commonly, 
the time intervals T between atoms are Poisson-distributed according to dV{T) = 
exp{—RT)RdT with an average time interval 1/R between them. Averaging the 
exponential in Eq. (|2.19| ) we get 

{piT))T = Sp , (2.20) 

where 

and N = R/'y is the dimensionless pumping rate already introduced. 

Implicit in the above consideration is the lack of knowledge of the actual value 
of the atomic state after the interaction. If we know that the state of the atom is 
s = ± after the interaction, then the average operator that transforms the cavity 
state is instead 

S{s) = {l + Lc/N)-'M{s) , (2.22) 

with M(s) given by Eq. (|]T|). 

Repeating the process for a sequence of k unobserved atoms we find that the 
initial probability distribution p becomes S'^p. In the general case this Markov 
process converges towards a statistical equilibrium state satisfying Sp = p, which 



has the solution [^, ^ for tt, > 1 

nf,m + Naq^ 



m=i (1 +nb)m + Nbqm 



The overall constant po is determined by I^^^LoPn — ^■ 
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3 Correlations 



After studying stationary single-time properties of the micromaser, such as the aver- 
age photon number in the cavity and the average excitation of the outgoing atoms, 
we now proceed to dynamical properties. Correlations between outgoing atoms are 
not only determined by the equilibrium distribution in the cavity but also by its ap- 
proach to this equilibrium. Short-time correlations, such as the correlation between 
two consecutive atoms are difficult to determine experimentally, because 

they require efficient observation of the states of atoms emerging from the cavity in 
rapid succession. We propose instead to study and measure long-time correlations, 
which do not impose the same strict experimental conditions. These correlations 
turn out to have a surprisingly rich structure (see Fig. |l]) and reflect global proper- 
ties of the photon distribution. In this section we introduce the concept of long-time 
correlations and present two ways of calculating them numerically. In the following 
sections we study the analytic properties of these correlations and elucidate their 
relation to the dynamical phase structure, especially those aspects that are poorly 
seen in the single-time observables or short-time correlations. 



3.1 Atomic beam observables 

Let us imagine that we know the state of all the atoms as they enter the cavity, 
for example that they are all excited, and that we are able to determine the state 
of each atom as it exits from the cavity. We shall assume that the initial beam 
is statistically stationary, described by the density matrix (^78|) , and that we have 
obtained an experimental record of the exit states of all the atoms after the cavity has 
reached statistical equilibrium with the beam. The effect of non-perfect measuring 
efficiency has been considered in several papers B3, B3| but we ignore that 



complication since it is a purely experimental problem. From this record we may 
estimate a number of quantities, for example the probability of finding the atom in 
a state s = ± after the interaction, where we choose + to represent the excited state 
and — the ground state. The probability may be expressed in the matrix form 

V{s) = u^'^M{s)p'' , (3.1) 



where M{s) is given by Eq. ( |2.12| ) and p° is the equilibrium distribution ( p.23| ). The 
quantity is a vector with all entries equal to 1, m° = 1, and represents the sum 
over all possible final states of the cavity. In Fig. ^ we have compared the behaviour 
of 'P(+) with some characteristic experiments. 

Since V{+) +V{—) = 1 it is sufficient to measure the average spin value 

{s)=V{+)-V{-) . (3.2) 



Since = 1 this quantity also determines the variance to be (s^) — (s)^ = 1 — (s)^. 
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Correspondingly, we may define the joint probability for observing the states of 
two atoms, Si followed S2, with k unobserved atoms between them. 



n(si,S2)=M°^5(s2)5^'5(si)/ , (3.3) 



where S and S{s) are defined in Eqs. ( |2.21| ) and ( |2.22| ). The joint probability 



of finding two consecutive excited outcoming atoms, Vo{+,+), was calculated in 
Il30|| . It is worth noticing that since S = S{+) + S{—) and Sp^ = we have 
J2si 'Pkisi, S2) = V{s2)- Since we also have u^^ L = u^^ {M — 1) = we find likewise 
that u^^ S = u^~^ so that J2s2'^kisi, S2) = V{si). Combining these relations we 
derive that 'Pfc(+,— ) = ,+), as expected. Due to these relations there is 

essentially only one two-point function, namely the "spin-spin" covariance function 

= Vk{+,+)+Vk{-,-)-Vk{+,-)-Vk{-,+) (3-4) 
= l-4n(+,-) . 

From this we derive the properly normalized correlation function 

- . (3-5) 

which satisfies — 1 < 7^ < 1. 

At large times, when k —>■ 00, the correlation function is in general expected to 
decay exponentially, and we define the atomic beam correlation length by the 
asymptotic behaviour for large k c:^ Rt 

Here we have scaled with R, the average number of atoms passing the cavity per 
unit of time, so that C,a is the typical length of time that the cavity remembers 
previous pumping events. 

3.2 Cavity observables 

In the context of the micromaser cavity, one relevant observable is the instantaneous 
number of photons n, from which we may form the average (n) and correlations in 
time. The quantum state of light in the cavity is often characterized by the Fano- 



Mandel quality factor which is related to the fluctuations of n through 
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^'-^^) ^ ■ ^^-^^ 

This quantity vanishes for coherent (Poisson) hght and is positive for classical light. 

In equilibrium there is a relation between the average photon occupation num- 
ber and the spin average in the atomic beam, which is trivial to derive from the 
equilibrium distribution 

(n) = = nb + NV{-) = rife + ~}'^^ , (3.8) 



where n is a diagonal matrix representing the quantum number n. A similar but 
more uncertain relation between the Mandel quality factor and fluctuations in the 



atomic beam may also be derived p9 |. 



The covariance between the values of the photon occupation number k atoms 
apart in equilibrium is easily seen to be given by 

(rm)fc = , (3.9) 



and again a normalized correlation function may be defined 



The cavity correlation length is defined by 



7f~expf--^) . (3.11) 



Since the same power of the matrix 5* is involved, both correlation lengths are 
determined by the same eigenvalue, and the two correlation lengths are therefore 
identical C,a = = ^ and we shall no longer distinguish between them. 



3.3 Monte Carlo determination of correlation lengths 

Since the statistical behaviour of the micromaser is a classical Markov process it is 
possible to simulate it by means of Monte Carlo methods using the cavity occupation 
number n as stochastic variable. 

A sequence of excited atoms is generated at Poisson-distributed times and are 
allowed to act on n according to the probabilities given by Eq. (PT^). In these 
simulations we have for simplicity chosen a = 1 and 6 = 0. After the interaction 
the cavity is allowed to decay during the waiting time until the next atom arrives. 
The action of this process on the cavity variable n is simulated by means of the 
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transition probabilities read off from the dissipative master equation ( p.l5| ) using a 



suitably small time step dt. The states of the atoms in the beam are determined 
by the pumping transitions and the atomic correlation function may be determined 
from this sequence of spin values by making suitable averages after the system has 
reached equilibrium. Finally the correlation lengths may be extracted numerically 
from the Monte Carlo data. 

This extraction is, however, limited by noise due to the finite sample size which 
in our simulation is 10^ atoms. In regions where the correlation length is large, it 
is fairly easy to extract it by fitting to the exponential decay, whereas it is more 
difficult in the regions where it is small (see Fig. This accounts for the differences 
between the exact numerical calculations and the Monte Carlo data in Fig. |l|. It 
is expected that real experiments will face the same type of problems in extracting 
the correlation lengths from real data. 

3.4 Numerical calculation of correlation lengths 

The micromaser equilibrium distribution is the solution of Sp = p, where S is the 
one-atom propagation matrix (|2.21|) , so that p^ is an eigenvector of 5* from the 
right with eigenvalue kq = 1. The corresponding eigenvector from the left is m° and 
normalization of probabilities is expressed as m° p° = 1. The general eigenvalue 
problem concerns solutions to Sp = up from the right and S = ku^ from the left. 
It is shown below that the eigenvalues are non-degenerate, which implies that there 
exists a spectral resolution of the form 

oo 

S = Y.Kep'u'~' , (3.12) 

with eigenvalues ki and eigenvectors p^ and from right and left respectively. The 
long-time behaviour of the correlation function is governed by the next-to-leading 
eigenvalue ki < 1, and we see that 

R^ = --^ . (3.13) 
log/ti 

The eigenvalues are determined by the characteristic equation detjS" — k} = 0, 
which may be solved numerically. This procedure is, however, not well-defined 
for the infinite-dimensional matrix S, and in order to evaluate the determinant we 
have truncated the matrix to a large and finite-size K x K with typical K ~ 100. 
The explicit form of S in Eq. (|2.21| ) is used, which reduces the problem to the 



calculation of the determinant for a Jacobi matrix. Such a matrix vanishes outside 
the main diagonal and the two subleading diagonals on each side. It is shown in 
Section |4.3| that the eigenvalues found from this equation are indeed non-degenerate, 
real, positive and less than unity. 
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The next-to-leading eigenvalue is shown in Fig. |1| and agrees very well with the 
Monte Carlo calculations. This figure shows a surprising amount of structure and 
part of the effort in the following will be to understand this structure in detail. 

It is possible to derive an exact sum rule for the reciprocal eigenvalues (see 
Appendix ^ , which yields the approximate expression: 

when the subdominant eigenvalues may be ignored. Here p„ is the equilibrium 
distribution Eq. ( |2.23| ) and -P„ = J2^2oPm is the cumulative probability. In Fig. | 
we compare the exact numerical calculation and the result of the sum rule, which is 
much less time-consuming to compute. 



4 Analytic preliminaries 

In order to tackle the task of determining the phase structure in the micromaser 
we need to develop some mathematical tools. The dynamics can be formulated in 
two different ways which are equivalent in the large flux limit. Both are related to 
Jacobi matrices describing the stochastic process. Many characteristic features of 
the correlation length are related to scaling properties for N ^ oo, and require a 
detailed analysis of the continuum limit. Here we introduce some of the concepts 
that are used in the main analysis in Section ^. 



4.1 Continuous master equation 

When the atoms have Poisson distributed arrival times it is possible to formulate 



the problem as a differential equation [41]. Each atom has the same probability Rdt 



of arriving in an inflnitesimal time interval dt. Provided the interaction with the 
cavity takes less time than this interval, i.e. r <^ dt, we may consider the transition 
to be instantaneous and write the transition matrix as Rdt{M — 1) so that we get 

^ = -^Lcp + RiM-l)p=-^Lp , (4.1) 



where L = Lq — N{M — 1). This equation obviously has the solution 

p{t) = e-^^'p . (4.2) 



Explicitly we have 
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[rib + l){n6n,m -{n + l)5„+i,m) + nb{{n + l)6n,m - n6n,m+i) 

(4.3) 



and 

1 dpn 



■{ub + l){npn - (n + l)Pn+i) - ribiin + l)p„ - np^-i) 



7 dt (4 4) 

-N{{aqn+i + hqn)pn - aqnPn-i - bqn+iPn+i) ■ 

The equilibrium distribution may be found by the same technique as before, 
writing the right-hand side of Eq. (|4.4|) as J„+i — Jn with 



Jn = {{nb + l)n + Nbqn)Pn - {nbn + Naqn)Pn-i , (4.5) 

and setting J„ = for all n. The equilibrium distribution is clearly given by the 
same expression (|2.23| ) as in the discrete case. 



4.2 Relation to the discrete case 

Even if the discrete and continuous formulation has the same equilibrium distribu- 
tion, there is a difference in the dynamical behaviour of the two cases. In the discrete 
case the basic propagation matrix is S^, where = (1 -|- Lc/N)~^M, whereas it is 
exp(— 7Lt) in the continuous case. For high pumping rate N we expect the two 
formalisms to coincide, when we identify k ~ Rt. For the long-time behaviour of 
the correlation functions this implies that the next-to-leading eigenvalues ki of S 
and Ai of L must be related by 1/C = 7-^1 — — -RlogKi. 

To prove this, let us compare the two eigenvalue problems. For the continuous 
case we have 

{Lc-N{M-l))p = Xp , (4.6) 

whereas in the discrete case we may rewrite 5*^ = up to become 

Lc-^{M~1)'^P = nQ--1^P . (4.7) 

Let a solution to the continuous case be p{N) with eigenvalue A(A^), making explicit 
the dependence on N. It is then obvious that p{N/k) is a solution to the discrete 
case with eigenvalue k determined by 
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A 



(4.8) 



As we shall see below, for A?" S> 1 the next-to-leading eigenvalue Ai stays finite or 
goes to zero, and hence Ki — ^ 1 at least as fast as 1/N. Using this result it follows 
that the correlation length is the same to 0{1/N) in the two formalisms. 

4.3 The eigenvalue problem 

The transition matrix L truncated to size {K -|- 1) x [K -|- 1) is a special kind of 
asymmetric Jacobi matrix 



r Aq + Bo -Bi 

-Ao + -B2 

-A A2 + B2 -Bs 



where 



-Ak-2 Ak-1 + Bk-1 —Bk 
-Ak-1 Ak + Bk 



(4.9) 



An = nb(n -Fl) Naqn+i , 



Bn = {nb+l)n + Nbqn 



(4.10) 



Notice that the sum over the elements in every column vanishes, except for the first 
and the last, for which the sums respectively take the values Bq and Ak- In our 
case we have Bq = 0, but Ak is non-zero. For Bq = it is easy to sec (using row 
manipulation) that the determinant becomes AqAi ■ ■ ■ Ak and obviously diverges 
in the limit of i^T — > 00. Hence the truncation is absolutely necessary. All the 
coefficients in the characteristic equation diverge, if we do not truncate. In order 
to secure that there is an eigenvalue A = 0, we shall force Ak — instead of the 
value given above. This means that the matrix is not just truncated but actually 
changed in the last diagonal element. Physically this secures that there is no external 
input to the process from cavity occupation numbers above K, a not unreasonable 
requirement. 

An eigenvector to the right satisfies the equation LkP — Ap, which takes the 
explicit form 

- ^n-lPn-l + {An + Bn)Pn " ^n+lPn+l = >^Pn ■ (4-11) 
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Since we may solve this equation successively for pi,p2, ■■■ ,Pk given po, it follows 
that all eigenvectors are non-degenerate. The characteristic polynomial obeys the 
recursive equation 



det(Lx - A) = {Ak + Bk - A) det(L;^_i - A) - Ak-iBk det{LK-2 - A) , (4.12) 



and this is also the characteristic equation for a symmetric Jacobi matrix with 
off-diagonal elements C„ = — Hence the eigenvalues are the same and 
therefore all real and, as we shall see below, non-negative. They may therefore be 
ordered = Aq < Ai < ■ ■ ■ < A/^. The equilibrium distribution ( |2.23| ) corresponds 
to A = and is given by 

P°=Pm%^=pg y"""^r forn = l,2,...,ir . (4.13) 

Notice that this expression does not involve the vanishing values Bq = = 0. 

Corresponding to each eigenvector p to the right there is an eigenvector u to the 
left, satisfying = Am"*", which in components reads 

AniUn - Un-l) + BniUn - Un+l) = \Un ■ (4-14) 

For A = we obviously have m° = 1 for all n and the scalar product ■ p^ = 1. 
The eigenvector to the left is trivially related to the eigenvector to the right via the 
equilibrium distribution 

Pn = PnUn ■ (4-15) 

The full set of eigenvectors to the left and to the right {u^,p^ | £ = 0, 1, 2, . . . , K} 
may now be chosen to be orthonormal ■ p^ = Si^fi, and is, of course, complete 
since the dimension K is finite. 

It is useful to express this formalism in terms of averages over the equilibrium 
distribution {fn)o = Yln=o fnPn- Then using Eq. ( [4.15| ) we have, for an eigenvector 
with A > 0, the relations 

{Un)o = , 

iOo = 1 , (4-16) 
{unu'^)o = for A 7^ A' . 
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Thus the eigenvectors with A > may be viewed as uncorrelated stochastic functions 
of n with zero mean and unit variance. 

Finally, we rewrite the eigenvalue equation to the right in the form of Ap„ = 
Jn - Jn+i with 



(4.17) 



Using the orthogonality we then find 



K 



A = XI Un{Jn - Jn+l) = {Bn{Un - Mn-l)^)c 



(4.18) 



ra=0 



which incidentally proves that all eigenvalues are non-negative. It is also evident 
that an eigenvalue is built up from the non-constant parts, i.e. the jumps of u„. 



4.4 Effective potential 

It is convenient to introduce an effective potential Vn, first discussed by Filipowicz 
et al. 1 27] in the continuum limit, by writing the equilibrium distribution ( 2.23|) in 
the form 



Pn 1 



(4.19) 



with 



1 ^ , Ubm + Naqr, 



m=l 



+ nb)m + Nbqn 



(4.20) 



for n > 1. The value of the potential for n = may be chosen arbitrarily, for 
example Vq = 0, because of the normalization constant 



n=0 



-NVn 



(4.21) 



It is, of course, completely equivalent to discuss the shape of the equilibrium distri- 
bution and the shape of the effective potential. Our definition of Vn differs from the 



one introduced in Refs. [27, 38 in the sense that our Vn is exact while the one in 



was derived from a Fokker-Planck equation in the continuum limit. 
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4.5 Semicontinuous formulation 



Another way of making analytical methods, such as the Fokker-Planck equation, 
easier to use is to rewrite the formalism (exactly) in terms of the scaled photon 
number variable x and the scaled time parameter 6', defined by |^ 



n 

^ ~ N ' (4.22) 
e = grVN . 

Notice that the variable x and not n is the natural variable when observing the 
field in the cavity by means of the atomic beam (see (|3.9| )). Defining Ax = 
and introducing the scaled probability distribution p{x) = Npn the conservation of 
probability takes the form 

CO 

p{x) = 1 , (4.23) 

x=0 

where the sum extends over all discrete values of x in the interval. Similarly the 
equilibrium distribution takes the form 

p\x) = i e-^^(^-) , (4.24) 

Z X 

with the effective potential given as an "integral" 

X 

V{x) = J2 D{x') , (4.25) 

with "integrand" 

D{ ) 1 ^''"^ ~'~ '^^^■^'^ (4 26) 
(1 + nb)x + hq{x) 

The transition probability function is q{x) = sin^ 9y/x and the normalization con- 
stant is given by 

= 4 = E e-^""^^) . (4.27) 



In order to reformulate the master equation ([4.4|) it is convenient to introduce 
the discrete derivatives A_|_/(x) = f{x + Ax)—f{x) and A_/(a;) = f{x)—f{x — Ax). 
Then we find 

ldp(x) A+ 

7 dt Ax-^^ ' ' ^ ^ 
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with 



J(x) = (x — (a — b)q(x))p(x) H (n^x + aq(x))— — p(x) . (4.29) 

N Ax 



For the general eigenvector we define p{x) = Np^ and write it as p{x) 
with u{x) = Un and find the equations 



Xp{x) 



Ax 



J{x) 



p^{x)u{x) 



(4.30) 



and 



J(x) = ^p°(a;)((l + n^)x + hq{x))^u{x) 
JM Ax 



(4.31) 



Equivalently the eigenvalue equation for u{x) becomes 



Xu{x) = {x-{a- b)q{x))^u{x) - 



[ribX + aq(x))— — u(x) 
Ax 



(4.32) 



As before we also have 



(u(x))c 



{u{xn 




1 



(4.33) 



where now the average over p^{x) is defined as (/(x))o = J2x Ax/(x)p°(x). As before 
we may also express the eigenvalue as an average 



A 1 N 7 / \N / A_U(X) 

A = ^ ( ((1 + nb)x + bq{x)) 



(4.34) 



Again it should be emphasized that all these formulas are exact rewritings of the 
previous ones, but this formulation permits easy transition to the continuum case, 
wherever applicable. 



4.6 Extrema of the continuous potential 

The quantity D{x) in Eq. ( |4.26| ) has a natural continuation to all real values of x as 
a smooth differentiable function. The condition for smoothness is that the change in 
the argument 6y/x between two neighbouring values, x and x + Ax is much smaller 
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than 1, or ^ ^fx. Hence for — > oo the function is smooth everywhere and 
the sum in Eq. ( |4.25|) may be replaced by an integral 

Vix) = r dx' D(x') , (4.35) 
Jo 

so that D{x) = V'{x). In Fig. |^ we illustrate the typical behaviour of the potential 
and the corresponding photon number distribution in the first critical region (see 
Section [5^ ). Notice that the photon-number distribution exhibits Schleich- Wheeler 



oscillations typical of a squeezed state 43 



The extrema of this potential are located at the solutions to q{x) = x; they may 
be parametrized in the form 

X = {a — b) sin^ , 

^ 1 * (".se) 



y/a — b I sin 01 

with < < oo. These formulas map out a multibranched function x{6) with 
critical points where the derivative 

j^u^) = Y'u^) = {a + nb{a-b)){q{x)-xq'{x)) 

^ ^ {{1 + nb)x + bq{x)){nbX + aq{x)) ' ' 

vanishes, which happens at the values of satisfying = tan0. This equation 
has an infinity of solutions, = 0^, k = 0, 1, . . ., with 0o = and to a good 
approximation 

for k = 1,2,..., and each of these branches is double-valued, with a sub-branch 
corresponding to a minimum {D' > 0) and another corresponding to a maximum 
{D' < 0). Since there are always k + 1 minima and k maxima, we denote the 
minima X2k{0) and the maxima X2k+i{0)- Thus the minima have even indices and 
the maxima have odd indices. They are given as a function of 6 through Eq. ([4.36|) 



when runs through certain intervals. Thus, for the minima of V{x), we have 

0fc < < (fc + l)7r, ek<6<oo, a-b> X2k{0) > 0, = 0, 1, . . . , (4.39) 

and for the maxima 

kiT < (f)< (j)k, oo > 6 > 6k, < X2k+iid) < a-b, k = 1, . . . (4.40) 
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Here 9k = <pk/\ sin (pkWo, — b is the value of 9 for which the fc'th branch comes into 
existence. Hence in the interval 9k < 9 < 9k+i there are exactly 2K + 1 branches, 
xq, xi, X2, ■ ■ ■ , X2K-1, X2K, forming the K + 1 minima and K maxima of V{x). For 
< 9 < 9o = l/\/a — b there are no extrema. 

This classification allows us to discuss the different parameter regimes that arise 
in the limit of — > oo. Each regime is separated from the others by singularities and 
are thus equivalent to the phases that arise in the thermodynamic limit of statistical 
mechanics. 



5 Phase structure 

We shall from now on limit the discussion to the case of initially completely excited 
atoms, a = 1, 6 = 0, which simplifies the following discussion considerably. 

The central issue in this paper is the phase structure of the correlation length as a 
function of the parameter 9. In the limit of infinite atomic pumping rate, N — > oo, 
the statistical system described by the master equation (|2.10|) has a number of 
different dynamical phases, separated from each other by singular boundaries in 
the space of parameters. We shall in this section investigate the character of the 
different phases, with special emphasis on the limiting behaviour of the correlation 
length. There turns out to be several qualitatively different phases within a range 
of 9 close to experimental values. First, the thermal phase and the transition to the 



maser phase at ^ = 1 has previously been discussed in terms of (n) |T|, ^ . The 
new transition to the critical phase at 9i ^ 4.603 is not revealed by (n) and the 
introduction of the correlation length as an observable is necessary to describe it. In 
the large flux limit (n) and ((An)^) are only sensitive to the probability distribution 
close to its global maximum. The correlation length depends crucially also on local 
maxima and the phase transition at 9i occurs when a new local maximum emerges. 
At ~ 6.3 there is a phase transition in (n) taking a discrete jump to a higher value. 
It happens when there are two competing global minima in the effective potential for 
different values of n. At the same point the correlation length reaches its maximum. 
In Fig. ^ we show the correlation length in the thermal and maser phases, and in 
Fig. the critical phases, for various values of the pumping rate A^. 

5.1 Empty cavity 

When there is no interaction, i.e. M = 1, or equivalently g„ = for all n, the be- 
haviour of the cavity is purely thermal, and then it is possible to find the eigenvalues 
explicitly. Let us in this case write 

Lc = {2nb + l)L3-{l + nb)L_-nhL+-^ , (5.1) 
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where 

(-^3)nm = ('^ + l) ^nm , 

These operators form a representation of the Lie algebra of SU(1,1) 

[L_,L+] = 2L3, [L3,L±] = ±L± . (5.3) 

It then follows that 

Lc = e^^+e-(^+"'')^- (L3 - i) e(i+"'')^-e-'-^+ , (5.4) 

where r = nb/{l + rii,). This proves that Lc has the same eigenvalue spectrum as 
the simple number operator L3 — |, i.e. A„ = n for n = 0, 1, . . .. Since M = 1 for 
r = this is a limiting case for the correlation lengths 7^^ = l/A„ = l/n for 6 = 0. 
From Eq. (^4.8|) we obtain k„ = 1/(1 + n/N) in the non- interacting case. Hence 
in the discrete case = — l/logK„ ^ N/n for N ^ n and this agrees with the 
values in Fig. |l] for n = 1, 2, 3 near r = 0. 

5.2 Thermal phase: < 6* < 1 

In this phase the natural variable is n, not x = n/N . The effective potential has no 
extremum for < n < 00, but is smallest for n = 0. Hence for ^ 00 it may be 
approximated by its leading linear term everywhere in this region 

iVV; = nlog^^^ . (5.5) 
rib + 9^ 

Notice that the slope vanishes for ^ = 1. The higher-order terms play no role as 
long as 1 — ^ 1/ V^, and we obtain a Planck distribution 

A^\-^{^X . (5.6, 



with photon number average 
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which (for 9 > 0) corresponds to an increased temperature. Thus the result of pump- 
ing the cavity with the atomic beam is simply to raise its effective temperature in 
this region. The mean occupation number (n) does not depend on the dimensionless 
pumping rate N (for sufficiently large A^). 
The variance is 

'^l = (-0 - {«)= = {">(1 + («)) = ''-'("I'y' . (5.8) 

and the first non-leading eigenvector is easily shown to be 

n — (n) 

Un= — , (5.9) 

which indeed has the form of a univariate variable. The corresponding eigenvalue is 
found from Eq. ( |4.34| ) Ai = 1 - or 

? = ^ . (5.10) 
Thus the correlation length diverges at ^ = 1 (for iV — > oo). 
5.3 First critical point: 9 = 1 

Around the critical point at 6' = 1 there is competition between the linear and 
quadratic terms in the expansion of the potential for small x 

V{x)=x\og^^^^^ + \x'-^ + 0{x') . (5.11) 
nb + 0^ 6 0^ + rib 

Expanding in — 1 we get 

V{x) = \^x + —^x' + 0{x\{e'-lf) . (5.12) 



Near the critical point, i.e. for (1 — (f'')\fN <^ 1, the quadratic term dominates, 
so the average value (x) as well as the width Ox becomes of (9(l/-\/iV) instead of 
0{\IN). 

Let us therefore introduce two scaling variables r and a through 



so that the probability distribution in terms of these variables becomes a Gaussian 
on the half-line, i.e. 
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-e 2 



(5.14) 



with 



(5,15) 



From this we obtain 



Ir) = a + = 



(5.16) 



For a = we have exphcitly 



12(1 + nb) 

TlN 



Q{nb + l) fl_l 

N \2 TT 



(5.17) 



This leads to the following equation for u{r) 

du d 

pu = rir — a)- — 

dr dr 



du 
dr 



{5.U 



where 



p = \^ 



' 3N 
l+rib 



' du 
dr 



(5.19) 



This eigenvalue problem has no simple solution. 

We know, however, that u{r) must change sign once, say at r = tq. In the 



neighbourhood of the sign change we have u ~ r — tq and, inserting this into ( |5.18 
we get ro = (a + v4 + f?)/2 and p = v4 + c? such that 



3N 

:i + nb)(4 + a2 



(5.20) 



5.4 Maser phase: 1 < ^ < ~ 4.603 

In the region above the transition at ^ = 1 the mean occupation number (n) grows 
proportionally with the pumping rate N, so in this region the cavity acts as a maser. 
There is a single minimum of the effective potential described by the branch Xq{6), 
defined by the region < < vr in Eq. (|4.36|) . We find for ^ 1 to a good 
approximation in the vicinity of the minimum a Gaussian behaviour 



23 



p\x) = ^^^^Me-f v"(-o)(.-.o)^ ^ (5.21) 

where 

Xo(l + nb) 

Hence for (^^ — 1)^/N ^ 1 we have a mean value (a;)o = xq and variance = 
l/NV"{xo). To find the next-to-leading eigenvalue in this case we introduce the 



scaling variable r = ^ NV"{xq){x — Xq), which has zero mean and unit variance for 
large N . Then Eq. ( |4.32| ) takes the form (in the continuum limit oo) 

Xu = {l-q'{x,)){r^-^] . (5.23) 

This is the differential equation for Hermite polynomials. The eigenvalues are = 
n(l — g'(xo)), n = 0, 1, . . ., and grow linearly with n. This may be observed in Fig. 
^. The correlation length becomes 

i = \-— = ] for < < TT . (5.24) 

1 — q [Xo) 1 — 0cot0 

As in the thermal phase, the correlation length is independent of N (for large N). 
5.5 Mean field calculation 

We shall now use a mean field method to get an expression for the correlation length 
in both the thermal and maser phases and in the critical region. We find from the 
time-dependent probability distribution ( [4.4[ ) the following exact equation for the 
average photon occupation number: 

^"^^""^ -Ar(g„+i)+n,-(n) , (5.25) 



or with Ax = 1/N 



7 dt 



1 d(x) 

^-^ = (q(x + Ax)) + UbAx - (x) . (5.26) 

7 dt 



We shall ignore the fiuctuations of x around its mean value and simply replace this 
by 
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ld(^ ^ q(^{^x) + Ax) + UbAx - (x) . (5.27) 
7 at 

This is certainly a good approximation in the hmit of iV — > oo for the maser phase 
because the relative fluctuation a^/ix) vanishes as 0(l/\/iV) here, but it is of du- 
bious validity in the thermal phase, where the relative fluctuations are independent 
of A^. Nevertheless, we find numerically that the mean field description is rather 
precise in the whole interval < ^ < ^i. 

The fixed point Xq of the above equation satisfies the mean field equation 

xo ^ q{xo + Ax) + UbAx , (5.28) 
which may be solved in parametric form as 

Xq — sin^ (p + UhAx , 

(/. (5.29) 

Jsin^ + (1 + rib) Ax 

We notice here that there is a maximum region of existence for any branch of the 
solution. The maximum is roughly given by — {k + 1)t:^N/{1 + n^). 
For small perturbations {x) — Xq -\- e we find the equation of motion 

-%^-{^-(i{.xo + Ax))e , (5.30) 
7 at 

from which we estimate the leading eigenvalue 

.^l-,'(..A.)^l--4f^l^. (5.31, 

Notice that A takes negative values in the unstable regions of 0. This eigenvalue 
does not vanish at the critical point 6 — 1 which corresponds to 

( ?,{l + nb) \^ 

TV J ' ^ 

but only reaches a small value 



A-2ji±p , (5.33) 
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which agrees exactly with the previously obtained result ( |5.20| ). Introducing the 
scaling variable a from (|5.13|) and defining i\} = (0/0o)^ "we easily get 



a = (V'=^-l)/^ , 

r = ^P , (5.34) 

and after eliminating ijj 



r 



2 



(5.35) 



p = a/4 + a2 , 
which agrees with the previously obtained results. 



5.6 The first critical phase: 4.603 ^ 9i < 9 < 02 ^ 7.790 

We now turn to the first phase in which the effective potential has two minima 
(xo, X2) and a maximum (xi) in between (see Fig. |^ in Section . In this case there 
is competition between the two minima separated by the barrier and for N —>■ 00 
this barrier makes the relaxation time to equilibrium exponentially long. Hence we 
expect Ai to be exponentially small for large N (see Fig. |^) 

= Ce"'?^ , (5.36) 



where C and rj are independent of A^. It is the extreme smallness of the subleading 
eigenvalue that allows us to calculate it with high precision. 

For large the probability distribution consists of two well-separated narrow 
maxima, each of which is approximately a Gaussian. We define the a priori proba- 
bilities for each of the peaks 

Po= E ^xp\x) = ^ , (5.37) 

0<x<xi ^ 

and 

P2= E Ax/(a;) = ^ . (5.38) 

X\ <x<oo 

The Z-factors are 
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Zo = EA^e-^^(^)^e-^^oJJ^ , (5.39) 



x=Q 

and 



£ Axe-^^(^)^e-^^^/-|^ , (5.40) 

x=xi V 2 

with ^ = ^0 + ^2. The probabihties satisfy of course Pq + P2 = 1 and we have 

p'^ix)=PopUx) + P2plix) , (5.41) 



where are individual probabihty distributions with maximum at xo,2- The over- 
lap error in these expressions vanishes rapidly for N — >■ 00, because the ratio P0/P2 
either converges towards or 00 for Vq 7^ V2. The transition from one peak being the 
highest to the other peak being the highest occurs when the two maxima coincide, 
i.e. at 6* ^ 7.22 at = 10, whereas for = cxd it happens at 6* ^ 6.66. At this 
point the correlation length is also maximal. 

Using this formalism, many quantities may be evaluated in the limit of large A^. 
Thus for example 

{x)o = Poxo + P2X2 , (5.42) 



and 

(^l = {{X - (x)o)')o = CTo'Po + + (Xo - X2fPoP2 ■ (5.43) 



Now there is no direct relation between the variance and the correlation length. 

Consider now the expression ( [4.3U| ), which shows that since Ai is exponentially 
small we have an essentially constant Jn, except near the maxima of the probability 
distribution, i.e. near the minima of the potential. Furthermore since the right 
eigenvector of Ai satisfies Y^xP{x) = 0, we have = J(0) = J(oo) so that 

< X < Xo , 
J{x) ^ <( Ji Xo < X < X2 , (5.44) 
X2 < X < oo . 

This expression is more accurate away from the minima of the potential, Xo and X2. 



Now it follows from Eq. (f4.31|) that the left eigenvector u{x) of Ai must be 
constant, except near the minimum xi of the probability distribution, where the 
derivative could be sizeable. So we conclude that u{x) is constant away from the 
maximum of the potential. Hence we must approximately have 
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u(x) ~ < (5.45) 

[ U2 Xi < X < OO . 

This expression is more accurate away from the maximum of the potentiaL 

We may now relate the values of J and u by summing Eq. ( [4.3U| ) from xi to 
infinity 

OO 

= J(xi) = Ai ^ Aa; p°(x)u(a;) ~ \1P2U2 . (5.46) 

X = Xl 



From Eq. ( 4.31 ) we get by summing over the interval between the minima 



N 1 ^"^ 1 
U2-u, = — ^ ^^^r-, ■ (5-47) 

The inverse probability distribution has for — * 00 a sharp maximum at the 
maximum of the potential. Let us define 



Then we find 



U2-uq = —— = — — . 5.49 



But uix) must be univariate, i.e. 



from which we get 



mqPo + U2P2 = , 
u^Po + UIP2 = 1 , 



\ Pq \ Zq 



U2 



(5.50) 



5 



(5.51) 





Inserting the above solution we may solve for Ai 

+ nbZo + . 
N Z0Z1Z2 
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or more explicitly 



v/^(^e-^(^-^«) + V^e-^(^-^^)) ■ (5.53) 



27r 

Finally we may read off the coefficients rj and C from Eq. (|5.36|) . We get 
and 

r^/y^ io^V.>V., (5,55) 
2n ^ [ for V2>Vo . 

This expression is nothing but the result of a barrier penetration of a classical 



statistical process |^ . We have derived it in detail in order to get all the coefficients 
right. 

It is interesting to check numerically how well Eq. ( ^.53|) actually describes the 
correlation length. The coefficient r] is given by Eq. ( |5.54] ), and we have numerically 
computed the highest barrier from the potential V{x) and compared it with an 
exact calculation in Fig. ^ The exponent t] is extracted by comparing two values 
of the correlation length, ^70 and ^90, for large values of (70 and 90), where the 
difference in the prefactor C should be unimportant. The agreement between the two 
calculations is excellent when we use the exact potential. As a comparison we also 
calculate the barrier height from the approximative potential in the Fokker-Planck 
equation derived in |2^, We find a substantial deviation from the exact value 
in that case. It is carefully explained in [^, ^ why the Fokker-Planck potential 
cannot be expected to give a quantitatively correct result for small n?,. The exact 
result (solid line) has some extra features at 6* = 1 and just below 6 = 61 c:^ 4.603, 
due to finite-size effects. 

When the first subleading eigenvalue goes exponentially to zero, or equivalently 
the correlation length grows exponentially, it becomes important to know the density 
of eigenvalues. If there is an accumulation of eigenvalues around 0, the long-time 
correlation cannot be determined by only the first subleading eigenvalue. It is quite 
easy to determine the density of eigenvalues simply by computing them numerically. 

In Fig. P we show the first seven subleading eigenvalues for = 50 and Ub = 
0.15. It is clear that at the first critical point after the maser phase {6 = 9i) there 
is only one eigenvalue going to zero. At the next critical phase {9 = 62) there 
is one more eigenvalue coming down, and so on. We find that there is only one 
exponentially small eigenvalue for each new minimum in the potential, and thus 
there is no accumulation of eigenvalues around 0. 
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6 Effects of velocity fluctuations 



The time it takes an atom to pass through the cavity is determined by a velocity 
filter in front of the cavity. This filter is not perfect and it is relevant to investigate 
what a spread in flight time implies for the statistics of the interaction between 
cavity and beam. To be specific, we consider the flight time as an independent 
stochastic variable. Again, it is more convenient to work with the rescaled variable 
6', and we denote the corresponding stochastic variable by i}. In order to get explicit 
analytic results we choose the following probability distribution for positive 

M«,/3) = :p^^^^'^e-'^" , (6.1) 
r(a + 1) 

with p = e/al and a = O'^/al - 1, so that {d) = 6 and {{^ - Of) = al Other 
choices are possible, but are not expected to change the overall qualitative picture. 
The discrete master equation ( |2.19|) for the equilibrium distribution can be averaged 
to yield 

(p(t + r))=e-^^-^(M(^))(p(t)) , (6.2) 

The factorization is due to the fact that p(t) only depends on {} for the preceding 
atoms, and that all atoms are statistically independent. The effect is simply to 
average q{'d) = sm'^{'dy/x) in M{{}), and we get 



1 - 1 + cos ^ arctan ' ^ ^ 



9^ \(yi \ 9 



(6.3) 



This averaged form of q{9), which depends on the two independent variables 
and 9^/x, enters in the analysis of the phases in exactly the same way as before. In 
the limit cr^ — we regain the original q{9), as we should. For very large ag and 
fixed 9, and (q) approaches zero. 

6.1 Revivals and prerevivals 

The phenomenon of quantum revival is an essential feature of the microlaser system 
(see |]15[-|jl9[, and 0]-|^). The revivals are characterized by the reappearance 



of strongly oscillating structures in the excitation probability of an outgoing atom 
which is given by Eq. (|3.1|): 



Vi+) = u'^M{+)p' = - qn-,m)p'n , (6.4) 

n 

where is the photon distribution ( p.23| ) in the cavity before the atom enters. Re- 
vivals occur when there is a resonance between the period in g„ and the discreteness 
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m n 



36[] . If the photon distribution in the cavity has a sharp peak at n = no with 
a position that does not change appreciably when 6 changes, as for example for a 
fixed Poisson distribution, then it is easy to see that the first revival becomes pro- 
nounced in the region of ^rev — 27r\/noN. For the equilibrium distribution without 
any spread in the velocities we do not expect any dramatic signature of revival, the 
reason being that the peaks in the equilibrium distribution p° (6*) move rapidly with 
6. In this context it is also natural to study the short-time correlation between two 
consecutive atoms, or the probability of finding two consecutive atoms in the excited 



level |30[. This quantity is given by 

Po(+,+) = M0'^M(+)(l + Lc/Ar)-iM(+)/ 



(6.5) 

n,m 



defined in Eq. ( |3.3| ). In Appendix |C| we give an analytic expression for the matrix 
elements of (1 -|- Lc/N)~^. In Fig. |10| we present V{+) and Vo{+,+) for typical 
values of N and rib. 

If we on the other hand smear out the equilibrium distribution sufficiently as a 
function of 9, revivals will again appear. The experimental situation we envisage 
is that the atoms are produced with a certain spread in their velocities. The sta- 
tistically averaged stationary photon distribution depends on the spread. After the 
passage through the cavity we measure both the excitation level and the speed of 
the atom. There is thus no averaging in the calculation of V{+) and 'Po(+, +), but 
these quantities now also depend on the actual value d for each atom. For definite- 
ness we select only those atoms that fall in a narrow range around the average value 
9, in effect putting in a sharp velocity filter after the interaction. The result for 
an averaged photon distribution is presented in Fig. ^ (lower graph), where clear 
signs of revival are found. We also observe that in 'Po(+;+) there are prerevivals, 
occurring for a value of 6 half as large as for the usual revivals. Its origin is obvious 
since in Vo{+,+) there are terms containing that vary with the double of the 
frequency of g„. 



6.2 Phase diagram 

The different phases discussed in Section |^ depend strongly on the structure of the 
effective potential. Averaging over 6 can easily change this structure and the phases. 
For instance, averaging with large ag would typically wash out some of the minima 
and lead to a different critical behaviour. We shall determine a two-dimensional 
phase diagram in the parameters 6 and ag by finding the lines where new minima 
occur and disappear. They are determined by the equations 
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The phase boundary between the thermal and the maser phase is determined by 
the effective potential for small x. The condition 6'^ = 1 is now simply replaced by 
("i?^) = 6"^ + ag = 1, which also follows from the explicit form of (q) in Eq. (|6.3| ). The 
transitions from the maser phase to the critical phases are determined numerically 
and presented in Fig. 0. The first line starting from 6' ~ 4.6 shows where the 
second minimum is about to form, but exactly on this line it is only an inflection 
point. At the point a about ag ~ 1.3 it disappears, which occurs when the second 
minimum fuses with the first minimum. From the cusp at point a there is a new 
line (dashed) showing where the first minimum becomes an inflection point. Above 
the cusp at point a there is only one minimum. Going along the line from point b to 
c we thus first have one minimum, then a second minimum emerges, and finally the 
first minimum disappears before we reach point c. Similar things happens at the 
other cusps, which represent the fusing points for other minima. Thus, solid lines 
show where a new minimum emerges for large n (~ A^) as 6 increases, while dashed 
lines show where a minimum disappears for small n as ae increases. We have also 
indicated (by dotted lines) the first-order maser transitions where the two dominant 
minima are equally deep. These are the lines where ^ and Qf have peaks and (n) 
makes a discontinuous jump. 

7 Finite-flux effects 

So far, we have mainly discussed characteristics of the large flux limit. These are the 
defining properties for the different phases in Section |^. The parameter that controls 
finite flux effects is the ratio between the period of oscillations in the potential and 
the size of the discrete steps in x. If g = sin^i^Oy/x) varies slowly over Ax = 1/A^, 
the continuum limit is usually a good approximation, while it can be very poor in 
the opposite case. In the discrete case there exist, for certain values of 9, states that 
cannot be pumped above a certain occupation number since g„ = for that level. 
This effect is not seen in the continuum approximation. These states are called 
trapping states |^ and we discuss them and their consequences in this section. 

The continuum approximation starts breaking down for small photon numbers 
when 6 > 2tt\/N, and is completely inappropriate when the discreteness is manifest 
for all photon numbers lower than A^, i.e. for 6 > 27tN. In that case our analysis 
in Section |^ breaks down and the system may occasionally, for certain values of 6, 
return to a non-critical phase. 
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7.1 Trapping states 



The equilibrium distribution in Eq. ( p. 231) has pecuhar properties whenever = 



for some value of m, in particular when Ub is small, and dramatically so when rib = 0. 



This phenomenon occurs when 6 = kTx^N/m and is called a trapping state. When 
it happens, we have p„ = for all n > m (for Ub = 0). The physics behind this 
can be found in Eq. ( [2.12 ), where M(— ) determines the pumping of the cavity by 



the atoms. If = the cavity cannot be pumped above m photons by emission 
from the passing atoms. For any non-zero value of nb there is still a possibility 
for thermal fluctuation above m photons and Pn 7^ even for n > m. The effect 
of trapping is lost in the continuum limit where the potential is approximated by 
Eq. ( [4.35| ). Some experimental consequences of trapping states were studied for 
very low temperature in and it was stated that in the range Ub = 0.1-1.0 no 
experimentally measurable effects were present. We, however, show below that there 
are clear signals of trapping states in the correlation length even for Ub = 1.0. 

7.2 Thermal cavity revivals 

Due to the trapping states, the cavity may revert to a statistical state, resembling 
the thermal state at 6' = 0, even if > 0. By thermal revival we mean that the 
state of the cavity returns to the 6 = thermal state for other values of 6. Even 
if the equilibrium state for non-zero 6 can resemble a thermal state, it does not at 
all mean that the dynamics at that value of 6 is similar to what it is at ^ = 0, 
since the deviations from equilibrium can have completely different properties. A 
straightforward measure of the deviation from the ^ = state is the distance in the 
norm 

Coo \ 1/2 

E[Pn(0)-Pn(^)f . (7.1) 

n=0 / 

In Fig. O we exhibit dL2[6) for = 10 and several values of nb. 



For small values of Ub we find cavity revivals at all multiples of v lOvr, which can 



be explained by the fact that sm{6yn/N) vanishes for n = 1 and iV = 10 at those 
points, i.e. the cavity is in a trapping state. That implies that pn vanishes for n > 1 
(for Ub = 0) and thus there are no photons in the cavity. For larger values of Ub the 
trapping is less efficient and the thermal revivals go away. 

Going to much larger values of 6 we can start to look for periodicities in the 
fluctuations in dL2{6). In Fig. [l^ (upper graph) we present the spectrum of periods 
occurring in di2{9) over the range < ^ < 1024. 

Standard revivals should occur with a periodicity of AO = 2ttsJ {n), which is 
typically between 15 and 20, but there are hardly any peaks at these values. On the 
other hand, for periodicities corresponding to trapping states, i.e. A^^ = nyJlO/n, 
there are very clear peaks, even though nb = 1.0, which is a relatively large value. 
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In order to see whether trapping states influence the correlation length we pre- 
sent in Fig. a similar spectral decomposition of ^ (6) (lower graph) and we find the 
same peaks. A more direct way of seeing the effect of trapping states is to study the 
correlation length for small n^. In Fig. |l^ we see some very pronounced peaks for 
small Tih which rapidly go away when rif, increases. They are located at 9 = irk \ 



N/n 

for every integer k and n. The effect is most dramatic when k is small. In Fig. 
there are conspicuous peaks at = 7rv^-{l/v^, 1/v^, 1, 2/^2}, agreemg 

well with the formula for trapping states. Notice how sensitive the correlation length 
is to the temperature when rib is small 150 



8 Conclusions 

We have thoroughly discussed various aspects of long-time correlations in the mi- 
cromaser. It is truly remarkable that this simple dynamical system can show such 
a rich structure of different phases. The two basic parameters in the theory are the 
time the atom spends in the cavity, r, and the ratio = -R/7 between the rate at 
which atoms arrive and the decay constant of the cavity. The natural observables 
are related to the statistics of the outgoing atom beam, the average excitation being 
the simplest one. We propose to use the long-time correlation length as a second 
observable describing different aspects of the photon statistics in the cavity. The 
phase structure we have investigated is defined in the limit of large flux, and can be 
summarized as follows: 

• Thermal phase, < 6* < 1. 

The mean number of photons (n) is low (finite in the limit N 00), and so is 
the variance cr^ and the correlation length ^. 

• Transition to maser phase, 6^ ~ 1. 

The maser is starting to get pumped up and S,-, (n), and cr„ grow like y/N. 

• Maser phase, 1 < 61 < 6*1 ~ 4.603. 

The maser is pumped up to (n) ~ A^, but fluctuations remain smaller, (T„ ~ 
V^, whereas ^ is finite. 

• First critical phase, 9i < 9 < 62 ^ 7.790. 

The correlation length increases exponentially with A^, but nothing particular 
happens with (n) and (T„ at 9i. 

• Second maser transition, 9 ^ 6.6 

As the correlation length reaches its maximum, (n) makes a discontinuous jump 
to a higher value, though in both phases it is of the order of A^. The fluctuations 
grow like A^ at this critical point. 

At higher values of 9 there are more maser transitions in (n), accompanied by 
critical growth of cr„, each time the photon distribution has two competing maxima. 
The correlation length remains exponentially large as a function of A^, as long as 
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there are several maxima, though the exponential factor depends on the details of 
the photon distribution. 

No quantum interference effects have been important in our analysis and the 
statistical aspects are purely classical. The reason is that we only study one atomic 
observable, the excitation level, which can take the values ±1. Making an analogy 
with a spin system, we can say that we only measure the spin along one direction. 
It would be very interesting to measure non-commuting variables, i.e. the spin in 
different directions or linear superpositions of an excited and decayed atom, and see 
how the phase transitions can be described in terms of such observables |]40| , |42 



Most effective descriptions of phase transitions in quantum field theory rely on clas- 
sical concepts, such as the free energy and the expectation value of some field, and 
do not describe coherent effects. Since linear superpositions of excited and decayed 
atoms can be injected into the cavity, it therefore seems to be possible to study co- 
herent phenomena in phase transitions both theoretically and experimentally, using 
the micromaser. 
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A Jaynes— Cummings with damping 



In most experimental situations the time the atom spends in the cavity is small 
compared to the average time between the atoms and the decay time of the cavity. 
Then it is a good approximation to neglect the damping term when calculating 
the transition probabilities from the cavity-atom interaction. In order to establish 
the range of validity of the approximation we shall now study the full interaction 
governed by the JC Hamiltonian in Eq. ( p.lD and the damping in Eq. (|2.14|) . The 
density matrix for the cavity and one atom can be written as 



p = p° (g) 1+ (g) (T^ + p+ ® Cr_ + p" (g) 0-_, 



(A.i; 



where p = ± ip^ and a± = (a^ ± iay)/2. We want to restrict the cavity part 
of the density matrix to be diagonal, at least the po part, which is the only part of 
importance for the following atoms, provided that the first one is left unobserved 
(see discussion in Section p7^) . Introducing the notation 



Pn 
Pn 



{n\po\n) , 
{n\pz\n) , 

(n|p_|_|n — 1) — {n — l|p_|n) 



the equations of motion can be written as 



(A.2) 



dt 

dpn 

dt 

dpn 

dt 



c „0 

m ; 



-^iVnpn + Vn + Ip^+i) - 7 II L^mP' 

i2g^{pl - pI_^ - Pn- Pn-l) -lY. LnmP\ 



± 

m ' 



(A.3) 



where 



c 



[{nb + l)n + nb{n + l)]6 

n,m [rib + l){n + 1) 

[nb{2n + 1) - i] 5n,m - {rib + + 1) 5n,m-i - nb^Jn{n - 1) • 

(A.4) 



It is thus consistent to study the particular form of the cavity density matrix, which 
has only one non-zero diagonal or subdiagonal for each component, even when damp- 
ing is included. Our strategy shall be to calculate the first-order correction in 7 in 
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the interaction picture, using the JC Hamiltonian as the free part. The JC part of 
Eq. ([A .31) can be drastically simplified using the variables 



Pi = Po - pr' - p: - pt' ■ (A.5) 

The equations of motion then take the form 



dpn ^ 7 
dt 2 



/ y IK-^nm ~ ^n—l,m—l J rm ' \-^nm ^n—l,m—l J rm 



dp 



± 



'^igV^Pn ~ X! [i^nm ~ -^n-l,m-l)Pm + i^nm + -^n-l,m-l)Pm )(A-6) 
— - - 2ig^pl--iY.Lnmpt ■ 

The initial conditions p^(0) = Pn-i, Pn(0) = ^Pn~i and p^(0) = are obtained from 

TT{p{O)\n){n\0Tl) = 2p^(0) = p„ , 
Tr(p(0)|n)(n|®^(ll-a,)) = p«(0) - p^(0) = , (A.7) 
Tr(p(0)|n)(n| = Ti {p{0)\n) {n\ ay) = . 



In the limit 7 — >■ it is easy to solve Eq. ( |A.6| ) and we get back the standard solution 
of the JC equations, which is 



Pait) = -Pn-icos{2gt^) , (A.8) 
p±(t) = -ipn^lsm{2gt^/n) . 



Equation ( [A.6|) is a matrix equation of the form p = [Cq — '~fCi)p. When Co and 



Ci commute the solution can be written as p{t) = exp(7Cit) exp(Cot)p(0), which 
is the expression used in Eq. ( p.l9| ). In our case Cq and Ci do not commute and 
we have to solve the equations perturbatively in 7. Let us write the solution as 
p{t) = exp(Cot)pi(t) since exp(Cot) can be calculated explicitly. The equation for 
pi{t) becomes 

^ = -ie-^°'C,e^^'pi{t) , (A.9) 
which to lowest order in 7 can be integrated as 
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pi(r) = -7 r dte-^°^Cie^"*p{0) + p{0) . (A.IO) 
Jo 



The explicit expression for exp(Cot) is 

/ 1 



cos{2gt^/n) ism{2gt^/n) \ , (A. 11) 

2 sm{2gt^/n) cos{2gty/n) 



and, therefore, exp(— Cot)Ci exp(Cot) is a bounded function of t. The elements of 
Ci are given by various combinations of L^^ and in Eq. ( |A.4| ) and they grow 
at most linearly with the photon number. Thus the integrand of Eq. (|A.1CI|) is of 
the order of (n) up to an nb-dependent factor. We conclude that the damping is 
negligible as long as 7r(n) -C 1, unless rib is very large. When the cavity is in a 
maser phase, (n) is of the same order of magnitude as iV = i?/7, so the condition 
becomes tR <^ 1. 



B Sum rule for the correlation lengths 

In this appendix we derive the sum rule quoted in Eq. (|3.14|) and use the notation 
of Section ^]3[ 

For Ak = the determinant det Lk becomes BqBi ■ ■ ■ Bk as may be easily 
derived by row manipulation. Since Ak only occurs linearly in the determinant 
it must obey the recursion relation det Lk = Bq - ■ ■ Bk + Ak det Lk-i- Repeated 
application of this relation leads to the expression 

K+l 

detLK=Y.Bo---Bk-iAk---AK . (B.l) 

fc=0 

This is valid for arbitrary values of Bq and Ak- Notice that here we define Bq - ■ ■ B^^i 
= 1 for k = and similarly A^. - ■ ■ Ak = 1 for k = K + 1. 

In the actual case we have Bq = Ak = 0, so that the determinant vanishes. The 
characteristic polynomial consequently takes the form 



det(L^ - A) = (-A) (Ai - A) ■ ■ ■ (A^^ - A) = -Di A + D2A' + ^(A^) , (B.2) 



where the last expression is valid for A ^ 0. The coefficients are 

D, = X,---Xk , (B.3) 



and 
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K 



Do 



1 



k=l 



(B.4) 



To calculate Di we note that it is the sum of the K subdeterminants along the 
diagonal. The subdeterminant obtained by removing the /c'th row and column takes 
the form 



-Ak^2 + Bk-i 



-Ak-i Ak + B 



K 



(B.5) 



which decomposes into the product of two smaller determinants which may be cal- 
culated using Eq. ( [B.l| ). Using that Bq = Ak = we get 



K 



Di = ^ Aq - ■ ■ ■ ■ ■ Bk 

k=0 



(B.6) 



Repeating this procedure for D2 which is a sum of all possible diagonal subdetermi- 
nants with two rows and columns removed {0 < k < I < K) , we find 



K-l K I 

-D2 = X! X! X! ^0' " ^k-lBk+l- ■ ■ Bm-lAm- ■ ■ Ai_iBij^i- ■ ■ Bk 
k=0 l=k+l m=k+l 



Finally, making use of Eq. ( 4.13 ) we find 
and 



Bi - ■ ■ Bk 

Po k=0 



PkPi 



R n K-l K I ^0^0 

rP 2^ 2^ 2^ r> 

yO k=0 l=k+l m=k+l ^rnPm 



Introducing the cumulative probability 



n-l 

n / y fm 5 
m=0 



(B.7) 



(B.8) 



(B.9) 



(B.IO) 
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and interchanging the sums, we get the correlation sum rule 



n=l ^" n=l -^nPn 

This sum rule is valid for finite K but diverges for ii' oo, because the equi- 
librium distribution p° approaches a thermal distribution for n ^ N. Hence the 
right-hand side diverges logarithmically in that limit. The left-hand side also di- 
verges logarithmically with the truncation size because we have A° = n for the 
untruncated thermal distribution. We do not know the thermal eigenvalues for the 
truncated case, but expect that they will be of the form A° = n + 0{n'^/K) since 
they should vanish for n = and become progressively worse as n approaches K. 
Such a correction leads to a finite correction to I]„l/A„. In fact, evaluating the 
right-hand side of Eq. ( B.ll ), we get for large K 



LtF-E ~ - - log(l + n^) . (B.12) 

n=l n=l "' n=l "' 



Subtracting the thermal case from Eq. (|B.11D we get in the limit oi K —>■ oo 



£U"A»j=„?A^« n j ■ (^■"' 

Here we have extended the summation to infinity under the assumption that for 
large n we have A„ ~ A°. The left-hand side can be approximated by ,^ — 1 in 
regions where the leading correlation length is much greater than the others. A 
comparison of the exact eigenvalue and the sum-rule prediction is made in Fig. |. 

C Damping matrix 

In this appendix we find an integral representation for the matrix elements of {x + 
Lc)^^, where Lc is given by Eq. ( |2.18| ). Let 

oo 

= E i^^nm + {Lc)nm)Wm , (C.l) 
r?i=0 

and introduce generating functionals v{z) and w{z) for complex z defined by 

oo oo 

v{z) = Y.Z^Vn , w{z) = Y.Z^Wn . (C.2) 

n=0 n=0 

By making use of 

oo 

^(^) = E + Lc)n7nZ^Wm , (C.3) 
n,m=0 
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one can derive a first-order differential equation for w{z), 

{x + n,{l-z))w{z) + {l + n,{l-z)){z-l)^^=v{z) , (C.4) 

which can be solved with the initial condition v{l) = 1, i.e. w{l) = 1/x. If we 
consider the monomial v{z) = VmZ™' and the corresponding w{z) = Wm{z), we find 
that 

Therefore {x+Lc)~^ is given by the coefficient of z"' in the series expansion ofwm{z). 
In particular, we obtain for nt = the result 

, r / m \ T(x + n)T(m - n + 1) . 

\ ^ J T{x + m + 1) 

where m > n. We then find that 



^o(+,+) = Ecos2(^?^v^ 

n=0 



"i)E 



m! A^r(A^ + n) 

; d'r(iv + m + 1; 



cos2((7rVm + , (C.7) 



where is the equilibrium distribution given by Eq. ( |2.23| ), and where x = N = 
R/j. Equation (|C.7|) can also be derived from the known solution of the master 
equation in Eq. ( |2.14| ) for rif, = |^^. For small and/or large x, Eq. ( |C.5| ) can 
be used to a find a series expansion in n^. 
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Figure 1: Comparison of theory (solid curve) and MC data (dots) for the correlation 
length (sample size 10^ atoms). The dotted and dashed curves correspond to subleading 
eigenvalues (k2,3 ) of the matrix S. The parameters are those of the experiment in Ref. J3l|/ . 
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Figure 2: Comparison of V{+) = 1—V{—) = l—{qn+i) with experimental data of Ref. [21] 



for various probability distributions. The Poisson distribution is defined in Eq. (2^), the 
thermal in Eq. ( |2.i(^ , and the micromaser equilibrium distribution in Eq. (2.2c). In the 
upper figure (N = R/'j = 1) the thermal distribution agrees well with the data and in 
the lower (N = 6) the Poisson distribution fits the data best. It is curious that the data 
systematically seem to deviate from the micromaser equilibrium distribution in the lower 
figure. 
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Figure 3: Monte Carlo data (with 10^ simulated atoms) for the correlation as a function 
of the separation k ^ Rt between the atoms in the beam for t = 25 fis (lower data points) 
and T = 50 fis (upper data points). In the latter case the exponential decay at large times 
is clearly visible, whereas it is hidden in the noise in the former. The parameters are those 
of the experiment described in Ref. ^Sl]. 
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Figure 4: Comparison of the sum in Eq. I \3.l4 ) over reciprocal eigenvalues (dotted curve) 
with numerically determined correlation length (solid curve) for the same parameters as in 
Fig. ||. The difference between the curves is entirely due to the suhdominant eigenvalues 
that have not been taken into account here. 
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Figure 5: Example of a potential with two minima xq,X2 and one maximum xi (upper 
graph). The rectangular curve represents the exact potential ( 4.2(\) , whereas the continu- 
ous curve is given by Eq. ( 4.230 with, the summation replaced by an integral. The value of 
the continuous potential at x = has been chosen such as to make the distance minimal 
between the two curves. In the lower graph the corresponding probability distribution is 
shown. 
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Figure 6: The correlation length in the thermal and maser phases as a function of 9 for 
various values of N. The dotted curves are the limiting value for N = oo. The correlation 
length grows as Vn near 9 = 1 and exponentially for 6 > ~ 4.603. 
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Figure 7: The logarithm of the correlation length as a function of 6 for various values 
of N (10, 20, ... , 100). We have nt, = 0.15 here. Notice that for 9 > Oi the logarithm 
of the correlation length grows linearly with N for large N. The vertical lines indicate 
eo = 1, 9i = 4.603, 02 = 7.790, 83 = 10.95 and 64 = 14.10. 
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Figure 8: Comparing the barrier height from the potential V{x) with the exact correlation 
length and the barrier from an approximate Fokker-Planck formula. 
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Figure 10: Upper graph: Probabilities of Ending one atom, or two consecutive ones, in the 
excited state. The Bux is given by N = 20 and the thermal occupation number isub = 0.15. 
The curves show no evidence for the resonant behaviour of revivals. Lower graph: Presence 
of revival resonances in equilibrium after averaging the photon distribution over 9. The 
same parameters as in Fig. [I^ are used but the variance in 9 is now given by ctq = 10. 
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Figure 11: Phase diagram in the 0-ae plane. The soUd hnes indicate where new minima 
in the effective potential emerge. In the lower left corner there is only one minimum at 
n = 0, this is the thermal phase. Outside that region there is always a minimum for 
non-zero n implying that the cavity acts as a maser. To the right of the solid line starting 
at ^ 2± 4.6, and for not too large ae, there are two or more minima and thus the correlation 
length grows exponentially with the flux. For increasing ag minima disappear across the 
dashed lines, starting with those at small n. The dotted hnes show where the two lowest 
minima are equally deep. 
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Figure 12: Distance between the initial probability distribution pn{0) andpn{d) measured 



bydL2{e) in Eq. (7.1). 
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Figure 13: Amplitudes of Fourier modes of di2{6) (upper graph) and ^{6) (lower graph) 
as functions of periods using N = 10, Ub = 1.0 and scanning < 9 < 1024. There are 
pronounced peaks at the values of trapping states: A9 = i^^fWjn. 
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